Hot electrons in a tunnel structure based on metal nanoclusters 
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We study the effect of temperature on the tunnel current in a structure based on gold clusters 
taking into consideration their discrete electronic spectra. We suggest that an overheating of electron 
subsystem leads to the disappearance of a current gap and gradual smoothing of current-voltage 
curves that is observed experimentally. 



The nanodispersed metallic systems are prospective 
object of nanotechnology Transport of electrical 

charge across a nanoscale tunnel junction is accompanied 
by many effects, such as the Coulomb blockade of the av- 
erage current, transfer of energy between electrons and 
ions, and consequent heating of the junction. In nanome- 
ter scale devices, electron transport can occur through 
well-resolved quantum states. If the temperature is in- 
creased, the Coulomb and quantum staircases of current 
are gradually smeared out by thermal fluctuations. 

Simple tunnel construction can be schematically rep- 
resented by the distinctive "sandwich" [T, 's']. It con- 
sists of a thick gold film (emitter) covered by a dielectric 
one (with dielectric constant e ~ 3). Disc-shaped [1] or 
spherical-like [3] gold clusters are self-organized on the 
dielectric layer. Also, a tip of STM is used in the capac- 
ity of the third electrode (collector). 

Some of the experimental features of the I — V curves 
were investigated in Ref. [4], however, the fact of smooth- 
ing of staircases for granule-molecule at low temperatures 
is still not understood. Such a smoothing is typical for 
molecular transistors Q. Moreover, the observed current 
gap decreases significantly as temperature increases from 
5 K to 300 K in structure based on disk-shaped cluster 
(Figure 2 in for disk 2R ^ 4nm). 

In this letter, the temperature features of the I — V 
curves are explained by overheating of electron subsys- 
tem. 

The number of atoms Nq {14, 10^} and {100, 600} 
correspond to gold discs of monatomic thickness and 
spheres whose radii vary in the range 2R {1, 8.5} and 
{1.4, 2.8} nm, respectively. For given cluster sizes, the 
condition L ^ R is fulfilled, where L is the free path 
length for the electrons in the bulk of a metal. 

The calculation of the electron spectrum in the cylin- 
drical and spherical wells of the mentioned sizes with fi- 
nite deepness yields different values for the spectrum dis- 
creteness in magic clusters Asp = e^^ — e^^ (see Figure 
1). In the nonmagic clusters the levels of lowest unoccu- 
pied states coincide with those of highest occupied ones, 
= ^HO at T = 0. 
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FIG. 1: Inset shows the specific difference between energies 
of lowest unoccupied electron state s^^ and highest occupied 
one in neutral discs (red) and spheres (blue) Aua^q 
T = 0. 



The characteristic Coulomb energy of charging is Eq = 
e^/C, where C is self-capacitance of single granule in a 
vacuum (in the case of a disc, the capacitance can be 
estimated as for the oblate spheroids of equal volume). 
The calculations of Ref. demonstrated that these C 
are too small for the width of the current gap to be ex- 
plained. The most obvious example is the case of a disc, 
since almost half of the disc surface contacts to the di- 
electric film. Therefore, for these granules we change 
C ^ (1 + e)C/2. Then, for discs and spheres we have 
Ec ^ {1.60, 0.21} and {1.82, 1.06} eV, respectively. We 
note that the value of the capacitance is sensitive to the 
shape of the granule surface, and even small deviation 
from the spherical shape can change significantly the ca- 
pacitance. 

The consequence of the phonon spectrum deformation 
of granules is the weakening of the electron-phonon inter- 
action within them: vp/R^ ojd, where is the electron 
velocity at the Fermi surface in the bulk, and cjd is the 
Debye frequency. This interaction can be so suppressed 
that the electron-electron interaction becomes the main 
mechanism for the dissipation of the energy, which is in- 



2 



emitter 



collector 




FIG. 2: Energy profile of structure for V > 0. /i^ is tiie elec- 
tron chemical potential of charged granule in external electric 
field in. 



jected to the particle. This additional energy results in 
the overheating of the electron subsystem, which is de- 
scribed by the Fermi statistics with some effective (en- 
hanced) temperature T^^, and the temperature of the ion 
subsystem only slightly changes Q . With the increase of 
the bias voltage V, the number of electrons, relaxing in 
the granule, increases significantly. 

Among them are all the electrons with energies in the 
interval erjV (see Figure 2) below the Fermi level of the 
granule {r]V is the fraction of the bias voltage on the 
granule) , since the "flow" of tunneling electrons increases 
from below lying levels, thereby, involving large number 
of conductivity electrons to the relaxation process. At 
the same time, channels of losses appear, which are re- 
lated to the generation of holes on the occupied levels 
and their subsequent recombination. The granule does 
not fragmentize at the significant overheating of the elec- 
tron subsystem, because the I — V curves are reproduced 
at the cyclic changes of the bias voltage 0, 0] . 

The estimate of the energy, which is pumped by the 
conductivity electrons to the granules of discontinuous 
films, is given in Ref. 01 (~ 0.2, 0.3 eV). This means 
that the experiments correspond to the Coulomb 

blockade regime in the region of current gap at the whole 
diapason of R and reasonable values of T^^. Also, the 
quantum ladder can be smeared out by the thermal fluc- 
tuations, 

Ec > Aep > fceT^'^'^ 

where Asp is the difference between discrete levels in the 
vicinity of the granule Fermi level, and Asp = Ae^ for 
magic clusters at T = 0. 

According to the simple model of Ref. [1] , we represent 
the emitter and the collector as the electron reservoirs 
with continuum spectrums, which are occupied in accor- 
dance with the Fermi-Dirac distribution with chemical 
electron potential fi^^ < and temperatures T®'^ equal 
to the thermostat one. In all cases the energy is counted 
off from the vacuum level. The electron chemical poten- 
tial /i^ of a granule in a quantum case can be defined by 
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FIG. 3: Calculated I — V curves of structure based on magic 
clusters: disc Au230 and sphere Au256- 1 - Au23o: dotted 
curve - T^'^'S = 5 K, solid (red) - T^'^ = 5 and T^^ = 2000 K. 
2 - Au23o: dotted curve - T^'^'^ = 300 K, solid (blue) - T^'^ = 
300 and T^^ = 2000 K. 3 - Au256: dotted curve - T^'"" = 30 K 



and = 2000 K, solid (green) - T^^^ = 300 and T 
2000 K. 



iff 



the normalization condition at a given temperature T^^: 



oo 

E 

p=l 



{1 + exp[(ep - AiS)/fcBTSj]} ^ = No, 



(1) 



The summation in ([T]) is performed for all one-particle 
states, A^o is the average number of conduction electrons 
in a granule. The spectrum of states is calculated in 
advance and, therefore, the chemical potential of neutral 
Auatq granules and its temperature dependence can be 
found from equation ([1]). 

The current flowing through a metallic quantuni_gran- 
ule (with restriction on its Coulomb instability [J]), is 
determined as P = or 



'''max 



Umax 

^P„{Z-Z}, (2) 



where Pn (the probability of the finding "in average" of n 
surplus electrons at the granule) is defined by the master 
equation in the stationary limit. In fact, the reduced cur- 
rent is calculated / = I/{ePoT^)^ where F®'^ are the tun- 
nel rates. In order to find values of P^^o/^o we use the 
recurrent relation. uJ^ are the partial electron "streams" 
from the last electrodes to the granule and in the oppo- 
site direction. For the comparison with the results of Ref. 
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[l|, 0] , the calculations are done for three temperatures of 
the collector and emitter = = 5, 30, 300 K, and 
also T^^ = T^ 2000 K. The values = = 1 and 
77 = 1/2 are used for all cases. 

Figure 3 shows calculated I—V curves for disc of radius 
R = 2nm (magic cluster AU230) and sphere of radius 
R = Inm (magic cluster Au256)- For low temperatures 
(feT^ff < A^f), the current gap width AVg = \Vo- \ + 
Vb+ is determined analytically by the conductance gap 
boundaries Vq- and Vo+. For example, Vq-^ is defined 
from the condition of absence of collector current of the 
direct I — V curve branche ( F > 0), and finally we have: 



AVa 



1 



1 



2 — 7] 1 + 77 



(3) 



where As = /i^ _ ^ho > q at T = and 77 is fixed as for 
y > 0. Calculated values of AVg are in a good agreement 
with the experimental values based both on spherical and 
disc-shape clusters. 

The calculation of the I — V curves and current gap 
can be done only numerically at A^B^elr — ^^f, when 
the larger part of the spectrum, compared to Asp, is 
responsible for the charge transfer. Our calculations show 
an evident dependence of / — V curves fiatness on the 
electron subsystem temperature. 

However, it order to obtain an agreement with ob- 
served I — V curves it is necessary to suggest that elec- 
trons in the emitter and collector are also heated up to 



some effective temperature, which is higher than the ther- 
mostat one. It is possible, because electrons (the current 
/ = IpA is provided by I/e ~ 10^ number of electrons 
per second) relax in generally on the free path length in 
the last electrodes. 

For the illustration, we present our result at Figure 3 
for the sphere at T^^ = 300 K and Tjj = 2000 K. Only 
by such a way, we can explain the fiattening of the I — V 
curves for the metallic cluster structures at low thermo- 
stat temperatures. With the increase of the bias voltage, 
the current fiow is accompanied by the increase of the 
electron gas temperature. 

In conclusion, we have calculated the I — V character- 
istics of structure based on magic clusters: disc AU230 
and sphere AU256 • We have suggested that the overheat- 
ing of electron subsystem leads to the disappearance of 
current gap and significant flattening of current-voltage 
curves. Our results are in agood qualitative agreement 
with experiment of Ref. [J, 
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